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Buried two dimensional electron gasses (2DEGs) have recently attracted 
considerable attention as a testing ground for both fundamental physics and 
quantum computation applications [IH3]. Such 2DEGs can be created by 
phosphorus delta (5) doping of silicon, a technique in which a dense and 
narrow dopant profile is buried beneath the Si surface. Phosphorous 5-doping 
is a particularly attractive platform for fabricating scalable spin quantum bit 
architectures [TH3] , compatible with current semiconductor technology. The 
band structure of the <5-layers that underpin these devices has been studied 
intensely using different theoretical methods |3HS]? but it has hitherto not been 
possible to directly compare these predictions with experimental data. Here 
we report the first measurement of the electronic band structure of a 5-doped 
layer below the Si(OOl) surface by angle resolved photoemission spectroscopy 
(ARPES). Our measurements confirm the layer to be metallic and give direct 
access to the Fermi level position. Surprisingly, the direct observation of the 
states is possible despite them being buried far below the surface. Using this 
experimental approach, buried states in a wide range of other material systems, 
including metallic oxide interfaces | |10H12| . could become accessible to direct 
spectroscopic investigations. 

In order to probe the electronic structure of the 5-layer states, samples were prepared 
by substitutional^ incorporating 1/4 of a monolayer (ML) of phosphorus atoms on a clean 
Si(OOl) surface and burying this structure under 2.0±0.6 nm of silicon [13]. Control samples 
were fabricated by similarly growing silicon on the clean surface, without the incorporation 
of a phosphorus layer. Fig. [T]a,e and b,f show ARPES measurements at two photon energies 
(36 and 113 eV) of a control sample and a 5-doped sample, respectively. Both samples show 
several similar structures that can be ascribed to silicon bulk and surface states, but the 
5-doped sample clearly shows sharp additional states that are missing on the control sample 
(enclosed in yellow). Most strikingly, we observe a state at the Fermi energy (Ej, < 0.2 eV), 
immediately confirming the predicted metallic nature of the 5-doped layers. Other additional 
states are found at a binding energy Ef, « 1 eV with a character similar to the bulk V valence 
band (VB) states. 

ARPES can probe the dimensionality of the observed states by acquiring data at different 
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FIG. 1: Band structure measurements of buried phosphorus 5-layers in Si(OOl). a and b 

ARPES measurements at hv = 36 eV for the control sample and 5-doped silicon respectively. The 
additional states due to the <5-layer are enclosed in yellow and marked, c and d; the photoemission 
intensity at normal emission, plotted as a function of k± for the same samples with the calculated 
bulk band structure [13] overlaid in red. <5-layer states are marked as before and surface states are 
enclosed in white in c. ARPES measurements taken at hv = 113 eV for the control sample e and 
the 5-doped sample f. The bulk Brillouin zone g and the definition of the axes k± , k x and k y , with 
the plane through the bulk X point extracted h. The green solid line along k x corresponds to the 
bulk directions seen in the ARPES measurements in a, b, e and f. 

photon energies [T3j . Varying the photon energy corresponds to a change in the wave vector 
perpendicular to the surface (k±) and while bulk states generally show dispersion as a 
function of k±, the two-dimensional states in the £-layer or genuine surface states do not. 
Fig. [I] c and d show the photoemission intensity in normal emission (k x = k y = 0) for the 
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control sample and £-doped sample, respectively, across the available photon energy range 
(14-200 eV). The horizontal axis has been converted from photon energy to kj_ assuming 
free electron final states [13] . The marked states in Fig. [I] b and f are clearly identified as 
two-dimensional by their lack of dispersion. Furthermore, they can be distinguished from 
the genuine surface states or surface resonances (emphasised by white contours at Ef> fs0.5 
and 2.0 eV) of Si(001), since these are present on both samples. Note that the surface 
states appear quite broad. This is ascribed to crystalline imperfection in the surface of the 
overlayer, grown at low temperature to minimise dopant segregation. 

The ARPES measurements in Fig. [T] c and d are plotted together with the calculated 
many-body bulk band structure for Si (after Ref. [14J) and the agreement is very good. The 
k± values for reaching the bulk X high symmetry point are found to correspond to photon 
energies of 36 and 113 eV. Although the bulk states are essentially unchanged by 5-doping, 
the introduction of a layer of n-type dopants has shifted the surface Fermi level toward the 
top of the band gap by ~ 0.23 eV. This results in a shift of the entire band structure (shallow 
bulk and surface states) to a higher binding energy as seen, for example, by comparing Figjl] 
a and b. 

The new electronic states in the 5-doped sample can be mapped in the plane parallel to 
the surface. Fig. [2] shows the dispersion in two high-symmetry directions (a and c), as well 
as the surface Fermi contour (b). The map directly reveals the periodicity of the dispersion 
parallel to the surface as corresponding to a (2 x 2) structure. This is to be expected because 
of the two rotational domains of (2 x 1) reconstruction present on a Si(001) surface. Note 
that this periodicity parallel to the surface does not necessarily imply the same periodicity 
in the <5-layer or in the bulk, it is anyway expected due to the well-known process of surface- 
umklapp scattering [13]. The Fermi surface contour appears to be streaked along the k x and 
k y directions. The streaks can be related to the different appearance of the state along the 
two cuts shown. In the k x = k y direction in Fig. [2] a, the 5-layer state appears as a sharp 
structure around k\\ = (k^ + ky) ' 5 = 0, whereas it is always observable along the k x = and 
k y = (Fig. [2] c) directions. The detailed reason for this is difficult to determine because 
the states are only barely visible at the Fermi energy and the detailed dispersion is not 
accessible. The dispersion could be anisotropic with a higher effective mass along k x and k y 
than along the diagonal direction. However, streaking would also be expected if step edges 
or the rotated 2x1 domains on the surface are restricted in the x and y directions. Similar 
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FIG. 2: Periodicity of 5-layer states, a and c ARPES measurements along k x = k y and k y = 0, 

respectively, for the 5-doped sample at hv = 113 eV with the <5-layer states enclosed in pink. The 
inset shows the detail of the <5-layer state with the calculated ir, 217 and 1A states (after Fig. 10 
d of Ref. [5]) overlaid, b A constant energy map at the Fermi level. The 2x2 surface unit cell is 
indicated by a green dashed box and the pale blue lines and arrows indicate where the band maps 
a and c have been extracted. The length of the surface reciprocal lattice vector g sur f is indicated. 

limited domain sizes in the 5-layer could also play a role but the present observations do not 
require the existence of new periodicities (such as 4 x n) in the delta layer [3]. 

Theoretical studies of 5-doped Si have primarily focused on the electronic states close to 
the Fermi energy, as these are the relevant states for transport jlHS]. A common prediction 
of the studies is the formation of impurity bands below the bulk conduction band minimum. 
The lowest lying states, called ir and 2T, derive from the projection of the conduction band 
minimum states onto the plane of the 5-layer. These states are non- degenerate, and different 
theoretical approaches have resulted in different calculated separation energies between them 
at k\\ = 0, referred to as valley splitting (see, for example, Refs. [HE]). Even for the present 
case of 1/4 ML, the Fermi level position of the IT band minimum and the valley splitting 
are dependent on factors such as the dopant distribution [5] . These values also vary between 
theoretical approaches, however for the present case of 1/4 ML binding energies of 300- 
500 meV for the ir band minima are typical 

Our observed metallic states can be assigned to the IF states and possibly also the 
2T states, depending on the size of the valley splitting. The observed occupation of the 
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FIG. 3: Band alignment and temperature dependence of the (5-layer states, a Low 

kinetic energy measurements of the Si 2p core level (after subtraction of a polynomial background) 
at normal emission for a range of photon energies, hv = 105.6 eV corresponds to a Ej, of just 
~1.3 eV. b Schematic representation of the band bending; the pale blue band indicates the silicon 
band gap. The red diamonds show the VBM position as extracted from the core level measurements 
(uncertainty in energy < 50 meV). The <5-layer state within the band gap is indicated in red 
(T=300 K) and blue (100 K). The (5-layer state at deeper binding energy is indicated schematically 
in yellow. The band bending for the same substrate, without the (5-layer is overlaid as a transparent 
white band, c Temperature dependent measurements of the valence band at normal emission and 
hv = 113 eV, and d the same data expanded to show the detail of the (5-layer state near the Fermi 
level and VBM. 

states is, however, much lower than predicted and appears to be < 200 meV. The most 
important reasons for this discrepancy is that the calculations assume an intrinsic Si bulk, 
in contrast to our strongly p-doped samples, and a symmetric potential around the 5-layer. 
It is nonetheless interesting to overlay calculated 5-layer states on the experimental data 
(inset of Fig. [2] a) - for the purpose of illustration, the IT, 2T and 1 A states adapted from 
Fig. 10 d of Ref. [3] are assumed. The Fermi level has been shifted by 140 meV in order to 
facilitate the comparison, but other than this the axes are the same. The width, depth and 
general shape of the measured state is consistent with the calculation. However, only one 
state is observed. This could be due to the limited experimental resolution or to a valley 
splitting that is slightly larger than calculated such that the 2r state remains occupied. 
In any event, it appears safe to assume that the next band, the 1A state, is sufficiently 
separated in energy to be unoccupied and hence to preclude its observation. 
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The 5- layer states at a binding energy of ~1 eV (Fig. [T] b) are spectroscopically quite 
similar to quantum-confined hole states [T5] but they are interpreted as valence band-derived 
impurity bands. A comparison to theoretical work is not possible at present for lack of 
published calculations in this energy range. These states might, however, give valuable 
information on the electronic structure of the 5-layers, even though they do not contribute 
to the actual transport. 

It is surprising that the electronic structure of the delta layers should be observable by 
ARPES at all. A simplistic view of photoemission suggests that the mean free path of the 
photoelectron is so short (<1 nm) that the photoemission intensity becomes vanishingly 
small for layers buried by just a couple of unit cells [IE]. Indeed, this is what happens 
to the photoemission intensity from the P 2p core level that is entirely attenuated after 6 
minutes of Si dosing, indicating negligible segregation of the dopants [HI [18] and allowing 
us to calibrate the Si overlayer thickness to be at least 1.4 nm for this evaporation time. 
The 5-layer states in the valence band and conduction band, on the other hand, are still 
observable for Si overlayer thicknesses in excess of 5 nm. The reason for this behaviour is the 
extended nature of these 5-layer states. While they are quite strongly confined to the region 
close to the 5-layer [7], their wave function still has a very small amplitude at the surface. 
This is sufficient to observe the states when they are resonantly enhanced by choosing the 
right photon energy. Such resonances are clearly observed in Fig. [T] d, where the states at 
the Fermi energy are most intense for photon energies that lead to emission near the bulk 
X points. This is consistent with the fact that the IT and 2r states are derived from the 
conduction band minimum. Such strong resonances in cross section are frequently observed 
for surface states [19] and surface 2DEGs [201 - 122] and permit the observation of states that 
have a wave function maximum far below the surface. 

Being able to directly probe the states of the 5-layer, relative to the surface and bulk states 
allows the band alignment in this region to be measured. As mentioned above, the position 
of the Fermi level at the surface is directly measured and the presence of the 5-layer leads 
to a shift of ~ 0.23 eV with respect to the surface of uniformly doped Si [23]. The position 
of the 5-layer states themselves where found to be independent from the thickness of the Si 
overlayer for a thickness range of nominally 2 to 5 nm (not shown). Note, however, that there 
is a significant uncertainty in the absolute Si overlayer thickness [13]. The band alignment 
on the bulk side of the 5-layer cannot be directly probed here for these thicknesses, but it 
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is reasonable to expect alignment back to the bulk Fermi level (which is straightforward to 
calculate [21] )• A schematic representation of the room temperature band alignment is given 
in Fig. [3]b. Note that this scenario is more complex than the situation usually assumed 
in calculations of the 5-layer's properties [1H9]. First of all, the bulk Si used is strongly 
p-doped, fixing the Fermi level on one side of the potential well. The Fermi level at the 
surface is largely dictated by the surface states (even though the presence of the 5-layer 
leads to a shift), such that the near surface region is seen to be n-type. Such details are 
important for the theoretical modelling of the ^-layer's electronic structure and ultimately 
for the operation of 5-layer based devices [IH3]. 

The schematic band bending near the surface can be confirmed by core level spectroscopy. 
To this end, Si 2p core level spectra were measured as a function of photon energy and thus 
electron kinetic energy E^ (Fig. [3] a) . Since the mean free path A of the photo-emitted 
electron has a dependence on E k [25], the sampling depth is also dependent on E k , with 
E k m 70 eV corresponding to a minimum in A of ~ 0.6 nm. A large increase in A can be 
achieved with very low [24J. It is possible to collect the Si 2p core level from E k < 2eV 
(sampling depth 2.5 nm [21]) to Ek > 70eV (sample depth m 0.6 nm), corresponding to 
depths within the region of strongest band bending, as illustrated in Fig. [3] b. Although 
there is significant uncertainty in estimating the sampling depth, this allows measured values 
of the band bending to be overlaid on Fig. |3]b (red diamonds). 

As an example of the possibilities opened by the direct measurement of £-layer states, 
we report the temperature dependence of the state close to the Fermi level. Fig. [3] c and d 
reveal that this state moves to higher binding energy as the temperature is reduced. Between 
100 and 300 K, the 5-layer state is seen to shift (reversibly) by around 80 meV, whilst the 
valence band maximum is essentially unchanged. This observed behaviour is somewhat 
different from the prediction of an almost constant position with respect to the Fermi level 
[H] , something that is probably caused by the aforementioned differences in bulk doping and 
shape of potential well between our experiment and most calculations. The important point 
is that the experimental approach presented here allows this temperature dependence to be 
directly measured, opening the possibility to test different calculational approaches. 

We emphasise that the experimental approach presented here is not limited to buried 
phosphorus dopants in silicon. 5-doping of other semiconductors (such as GaAs [26J and 
Ge [27]), boron (p— type) 5-layers [25] as well as metal oxide interfaces [TQ1 [TTj, [T6] could be 
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similarly studied. In addition to a plethora of such materials, a wealth of further information 
should also be extractable from such experimental data; we have already been able to discuss 
the temperature dependence to the 5-layer pinning, but more detailed studies should reveal 
parameters such as the effective mass [6] and the strength of the electron-phonon interaction 

[221 ED]. 
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Supplementary Information 
Experimental details 

ARPES measurements were performed at beamline 14 at MAXlab, using a Specs Phoibos 
electron analyser equipped with an imaging detector. The detector's angular imaging direc- 
tion corresponds to k x . Mapping orthogonal to this (k y ) is achieved by rotating the sample. 
The beamline has photon energy range 14 to 200 eV, but the extremely high cross-section 
of the Si 2p core levels (Eb ~ 100 eV) results in core-level excitation from second order light 
hindering measurements at photon energies of ?^100 eV. The beamline and analyser were 
set such that the energy and angle resolutions were ^50 meV and ^0.3°, respectively. The 
binding energies quoted are relative to a metal Fermi level. All measurements were carried 
out at room temperature, except those in the temperature scan of Fig. 3. For measuring 
the Si 2p core level at very low kinetic energy, as reported in Fig. 3 a, it was necessary 
to apply a bias voltage to the sample. Comparing the intensity of data taken at different 
photon energies (Fig. 1 c and d) necessitates a normalization of the spectra. To present the 
data in Fig. 1, each normal emission spectrum was normalized to the the total intensity of 
the detector image. 

Two types of samples are prepared, the 5-doped samples and the control samples. The 
starting point for both is a p-type Si(001) (0.01 Q cm) wafer. This is prepared in vacuo 
by degassing at 650°C followed by annealing to 1200°C, producing a clean and well-ordered 
surface with two rotational 2x1 domains. In order to prepare the control sample, an over- 
layer (psI.4 to 5 nm) of Si is then evaporated at a rate of 0.33 nm/min onto the surface 
held at 250°C. The deposited overlayer is not sufficiently well-ordered for the band structure 
to be observed with ARPES. However, a subsequent anneal to ^450° C reorders the over- 
layer such that bulk-like bands become visible. The surface states appear suppressed and 
broad compared with the initial surface (prior to overlayer growth) but the surface and bulk 
electronic structure of the samples prepared in this way is otherwise similar to those of the 
initial Si(001)-(2 x 1) surface. 

The 5-doped samples are prepared by depositing phosphorus on the clean Si(001) surface 
prior to the growth of a Si layer. To this end, phosphine gas, PH 3 , is dosed (1.4 L) with 
the substrate at room temperature. The sample is then annealed to 350°C for 60 s to 
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substitutionary incorporate phosphorus, P, into the surface layer. This procedure leads to 
a coverage of 1/4 monolayer of phosphorus atoms [S1,S2]. A Si layer is then grown on top 
and annealed as before, creating the buried 5-layer. 

The thickness of the final Si layer is estimated by monitoring the intensity of the P 2p 
core level, measured at a photon energy of 170 eV, as a function of Si dose. The com- 
plete suppression of intensity from this level is estimated to correspond to a minimum layer 
thickness of 1.4 nm. While the absolute thickness is error-prone, the relative thickness of 
the preparation is well controlled, and samples with overlayers of 150% and 250% of this 
minimum thickness were prepared. The 5-layer states are observed for all of these samples 
and their position is unchanged. 

Probing bulk states using ARPES 

For the reader not familiar with ARPES, we briefly summarise some of the main relevant 
points for photoemission from bulk states. The photoemission process for bulk electrons 
changes the perpendicular component of the electron's wave vector, k±, upon crossing the 
surface barrier. Indeed, strictly spoken k± is not even well-defined near the surface. Without 
the knowledge of the final state dispersion in the solid, the initial k± of the photoelectrons 
cannot be determined. 

A frequently used approach to solve this problem is the assumption of free electron final 
states [S3]. To consider a simplified case, we can assume that the effective mass of the final 
state band in the crystal is equal to the electron mass (m), and consider only photoemission 
along the surface normal (i.e. k x — k y — 0). In this case, it is straightforward to show that 
k± PS \J (2m/h 2 )[V(j + hv — Ef, — $] , where Vq is the inner potential, hv is the photon energy, 
Ej, is the binding energy of the electron and $ is the workfunction. All of these parameters 
are known [S4] and this approach has been used to convert the horizontal axis in Fig. 1 c 
and d from photon energy to k±. This conservation also identifies the photon energies of 36 
and 113 eV to correspond to emission from the bulk X states at normal emission. 

The wave vector parallel to the surface k\\, on the other hand, is conserved in the pho- 
toemission process such that k\\ in the panels of Fig. 1 a, b, e, f and Fig. 2 a and c actually 
correspond to the directions indicated in the figures. For surface states, surface resonances 
and two-dimensional 5-layer states, k\\ is the only relevant quantum number but for the 
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Supplementary Figure 1: Bulk Brillouin zones of Si(OOl). a Silicon Brillouin Zone (BZ) 
defining the axes k x , k y and k±. The bulk Brillouin zone is shown, with the axes definitions 
corresponding to the principle vectors of the primitive surface unit cell. Two planes with fixed k± 
corresponding to the BZ zone boundary and centre respectively are also shown, b several of the 
neighbouring BZs are also included such that the continuation of these plans into the 
neighbouring zones can be seen. The surface (2 x 2) reciprocal unit cell is overlaid in green 
(dashed square) , such that the implications of momentum exchange due to surface umklapp 
processes are apparent. 



photoemission from bulk states it is important to note that k± changes along the scan. 
Consider Supplementary Fig. la: For a photon energy of 36 eV or 113 eV, the plane probed 
by ARPES is curved, with the bulk X-point corresponding to normal emission. The bulk 
direction probed in Fig. 1 a, b, e, f and 2 c can be approximated to the bulk X — K direction 
and Fig. 2 a to the X — W direction. It is important to notice that when we refer to T, X, 
K and W, we refer to the high symmetry points in the bulk Brillouin Zone, as indicated in 
Supplementary Fig. 1 - this should not be confused with the 2-dimensional Brillouin Zone 
indices, for example in Ref. [S5]. 

Another important point is that, due to the two-dimensional periodicity, k\\ is only denned 
plus or minus a surface reciprocal lattice vector g sur f (and we choose to define the x and 
y axes to be in the 110 and 110 directions, such that g sur f is in the directions ±x and 
±y). To see the consequences of this for the photoemission of bulk states from a Si(001) 
surface, consider Supplementary Fig. lb which shows an extended Brillouin zone scheme. 
Variable photon-energy photoemission at normal emission probes the T — X direction of 
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the bulk band structure and using the free electron final state model, the positions of the 
high symmetry points are indicated in Fig. 1 c and d. When the bulk X point is reached, 
however, adding a surface reciprocal lattice vector (within the red plane) folds in bulk states 
from the T point of the neighbouring zone. This so-called surface umklapp process is the 
reason why bulk T states are frequently observed at photon energies corresponding to the 
bulk X points and vice versa [S6]. This is also evident in the present data, e.g. in Fig. 2 
c for the bulk X point reached at k± m 3.5 A" 1 where a replica of the bulk valence band 
maximum is clearly observed. 

Finally, we note that the surface umklapp processes in this particular case are not 
limited to the (1 x 1) reciprocal lattice of Si(OOl) but to the two rotated domains of a 
(2x1) reconstruction. This alone would give rise to the observed (2 x 2) periodicity of the 
5-layer states in Fig. 2 b. This observed periodicity does therefore not necessarily reflect 
the periodicity in the phosphorus 5-layer. 
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